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Quantum Information processing by NMR with small number of qubits is well established. 
Scaling to higher number of qubits is hindered by two major requirements (i) mutual coupling 
among qubits and (ii) qubit addressability. It has been demonstrated that mutual coupling can be 
increased by using residual dipolar couplings among spins by orienting the spin system in a liquid 
crystalline matrix. In such a case, the heteronuclear spins are weakly coupled but the homonuclear 
spins become strongly coupled. In such circumstances, the strongly coupled spins can no longer 
be treated as qubits. However, it has been demonstrated elsewhere, that the 2^ energy levels of a 
strongly coupled N spin-1/2 system can be treated as an N-qubit system. For this purpose the various 
transitions have to be identified to well defined energy levels. This paper consists of two parts. In 
the first part, the energy level diagram of a heteronuclear 5-spin system is obtained by using a newly 
developed heteronuclear z-cosy (HET-Z-COSY) experiment. In the second part, implementation of 
logic gates, preparation of pseudopure states, creation of entanglement and entanglement transfer 
is demonstrated, validating the use of such systems for quantum information processing. 

I. INTRODUCTION 



Theoretically, future quantum computers can simulate physical systems and solve certain problems 
more efficiently than classical computers [1-5]. This possibility has excited several researchers to 
strive for experimental realization of quantum computers. Apart from trapped ions and cavity-QED 
systems, Nuclear Magnetic Resonance (NMR) has proved as a suitable physical system for quantum 
information processing (QIP) [6,7]. In fact, NMR has emerged as the most successful of all systems, 
demonstrating several quantum algorithms including Shor's factorization algorithm [8-14]. However, 
scalability is one of the hurdles that NMR-QIP faces today. 

Originally it was conceived that liquid-state NMR will use as qubits, spin-1/2 nuclei having different 
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Larmor frequencies and weakly coupled to each other by indirect J-couplings [15-17]. Since J- 
couplings mediate through covalent bonds, they become vanishingly small beyond 4-5 bonds. Such 
systems are thus difficult to scale beyond 7-8 qubits. Attempts are going on to explore other couplings 
such as the direct dipolar couplings between spin- 1/2 nuclei, which are larger in magnitude, depend 
on the inter-nuclear distances, and not on the number of intermediate covalent bonds. However, in 
liquids, molecules undergo fast isotropic reorientations, thereby averaging out all dipolar interactions 
( homonuclear, heteronuclear, intermolecular and intra molecular). In solids, the molecules are static, 
and all the dipolar interactions are retained yielding broad spectra, rendering such systems unsuitable 
for quantum information processing. However, when molecules are dissolved in liquid crystal solvents, 
they have sufficient translational motion but get partially oriented and have anisotropic rotational 
motion. In such cases, while the intermolecular dipolar interactions are averaged to zero, the intra- 
molecular dipolar interactions survive which are scaled down from several KHz to a few KHz. Such 
systems yield well resolved spectra with finite number of sharp lines [18]. In such cases the spins 
have considerable dipolar couplings and it has been demonstrated that they can be used for quantum 
information processing [19-22]. 

However, in cases where the differences in Larmor frequencies or chemical shifts are not large 
compared to the dipolar couplings, the spins become strongly coupled. The eigenstates of the 
system are not products of the eigen-statcs of various spins but linear combinations of the product 
states, and spins can no longer be identified as a qubits [23]. Still, it has been demonstrated that 
the 2^ energy levels of the N spin- 1/2 system can be treated as 2^ states of a N-qubit system, and, 
can be used for quantum information processing [24]. Quantum Information processing has already 
been demonstrated in homonuclear 2,3 and 4-spin systems oriented in liquid crystal matrices [24]. 
In the homonuclear systems, it was found that if the spins are in symmetric chemical environment, 
the eigenstates break up into symmetry groups, leading to reduced number of qubits and other 
complications with degenerate transitions from different symmetry manifolds [21,22]. It is therefore 
desirable that the homonuclear spins are in non-symmetric chemical environment. Furthermore, it 
is advantageous to include heteronuclear spins in the system. Here, we extend our earlier work by 
including a heteronuclear spin with four homonuclear spins in non-equivalent chemical environment, 
all dipolar coupled to each other. We use the molecule l-chloro,2-flouro benzene (dissolved in liquid 
crystal solvent ZLI-1132), which yields a 5-qubit system in which the four protons are strongly 
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coupled to each other but weakly coupled to the fluorine. The Hamiltonian for the dipolar interaction 
among the strongly coupled protons (li) is, 

Hn = 27: Y: i^:;"(3/,,/,,-Ii.Ij), (1) 

i,j{i<j) 

where D°"=j^^ij is the scaled dipolar interaction, scaled by the order parameter [18]. The 
Hamiltonian for dipolar interaction between the protons and the fluorine(S) is 

Hd = 27r J2D^sh^S,, (2) 

i 

The full Hamiltonian of the system including Zeeman, dipolar and J-coupling terms, is given by: 

i 

+27r D^p{3IiJj,-iaj) + 27rYD''i:hIziS, + 27r ^ J^jhlj + 27rY J hshiS,. (3) 

The equilibrium spectrum of the protons and the fluorine is given in Fig. 1. In order to use 
this spin system for QIP, one needs to identify the energy levels of the system and to assign each 
transition to a pair of these levels. There are two possible methods. One well-known method is 
to numerically diagonalize the Hamiltonian of the system, with some initially guessed parameters 
(chemical shifts, dipolar and J-couplings) and iteratively fit the calculated and observed spectra 
[18]. The second method is the use of Z-COSY two-dimensional (2D) experiment which identifies 
connected transitions and draws automatically the energy-level diagram [25]. So far Z-COSY has 
been used only for homonuclear spins. This paper consists of two parts. In the first part (section 
II) we introduce a heteronuclear Z-COSY experiment and utilize it to obtain the complete energy 
level diagram of this dipolar coupled 5 spin system. In the second part (section HI) , we demonstrate 
that this spin system can be used as a 5-qubit system by preparing pseudopure states, performing 
C"-NOT operations and controlled SWAP operations, and creating and transferring entanglement 
using transition selective pulses. 

II. HET-Z-COSY 

The pulse sequence of a Z-COSY experiment is 90 — ti — a — r — /9 — ^2, where only the longitudinal 
magnetization is retained during the interval r by either phase cycling or by a gradient pulse [25,26]. 
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The 90'' pulse converts the equihbrium z-magnetization into coherences, which are frequency labeled 
during the period ti. The small angle a pulse ensures that each cross-section parallel to CU2 from 
the resulting 2D spectrum is equivalent to a one-dimensional (ID) experiment in which the peak 
corresponding to the diagonal is selectively inverted. The directly connected transitions to the 
inverted transition are finally measured in the linear regime by a small angle f3 pulse. The Z-COSY 
experiment has been successfully used to construct the energy level diagram in oriented homonuclear 
spin systems [25,27]. Here the Z-COSY experiment is extended to include a heteronucleus and we 
call this experiment as HET-Z-COSY. 

For "n" heteronuclear systems HET-Z-COSY experiment requires "n" 2D experiments where each 
experiment observes one of the heteronucleus during ^2- The pulse sequence for HET-Z-COSY 
experiment is 90 — — a — {Gz + kti) — P — t2, where the 90° and a pulses are applied on all 
the nuclei, and f3 is applied only on the observed nucleus. In the present case, two experiments 
have to be performed for detecting protons and fluorine respectively. Gz is the gradient pulse which 
retains only the the longitudinal magnetization. Besides creating multi-spin zz terms [25], the a- 
pulse creates all coherences including zero-quantum coherences. The gradient pulse destroys all 
other coherences as desired, except the homonuclear zero-quantum coherences. These zero-quantum 
coherences get partially converted into single quantum coherences by the last P pulse causing the well 
known zero-quantum interference in the 2D spectroscopy. To avoid this we have used the method 
of zero-quantum shifting [26] by incrementing the r delay as (Gz-|-kti), where k=2 was used. This 
shifts the zero-quantum peaks in the cui dimension which are eventually removed by symmetrizing 
the spectrum [25] . 

The 2D spectrum of HET-Z-COSY on the systems of oriented 1-chloro 2-fiouro benzene is shown 
in Fig. 2. The experimental spectrum when protons are observed in t2, is given on the left side 
((a) and (c)), while the experimental spectrum when fluorine is observed is given on the right side 
((b) and (d)). The spectrum in flgure 2(a) has the diagonal and connectivity information between 
the proton transitions, while figure 2(c) has the connectivity information of the fluorine transitions 
to the proton transitions. Figure 2(b) has the diagonal and connectivity information between the 
fluorine transitions, and flgure 2(d) has the connectivity information of the proton transitions to 
the fluorine transitions. Every connectivity is mutual such that if transition 'a' is progressively 
connected to transition '6', then transition "b' is also progressively connected to transition 'a'. Hence 
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the cross peaks should be symmetric about the diagonal [25]. The symmetrized spectrum (free of 
zero-quantum and other artifacts) is given in Fig. 3. A connectivity matrix is constructed from the 
symmetric 2D spectrum, by taking cross-sections parallel to lo^ at each transition along lo\. Prom this 
connectivity matrix the complete energy- level diagram can be constructed following the method given 
in reference 32 [25,27]. The complete energy level diagram constructed by the above procedure is 
given in fig. 4. All the steps including construction of the juxtaposed 2D spectrum, symmetrization, 
construction of the connectivity matrix and the construction of the complete energy-level diagram 
are automated by a program written in Matlab [28] . 

In this 5-spin system, there are 2^ = 32 energy levels. The proton transitions connect energy 
levels within two separate domains (A and B) of 16 energy levels each, corresponding respectively 
to -1-1/2 and -1/2 states of the fluorine spin (figure 4). The fluorine transitions connect the levels 
from A domain to the B domain. The various transitions observed in the spectrum are given serial 
numbers (1-60 for proton and 61-82 for the fluorine transitions). Except the four extreme eigenstates 
all others are linear combinations of the product states of the 4 protons. The extreme eigenstates 
are lA, 16A, IB and 16B respectively corresponding to the eigenstates \aaaaoL)^ \a(5(5(5(5), \(5aaaa) 
and where the first state corresponds to the fiuorine spin and last four to the protons. 

III. QUANTUM INFORMATION PROCESSING USING TRANSITION SELECTIVE PULSES 

Fig. 5 shows the labehng of the various eigenstates of the 5-spin system to 5-qubit labels. While 
labels can be chosen for optimum execution of various logical operations, here wc have retained a 
conventional labeling [29]. Transition selective (7r)-pulse has a narrow bandwidth and when tuned 
to the resonance frequency of a specific transition in the spectrum, it exchanges the amplitudes 
of the two eigen-states which are connected by that transition [23] . Such pulses have already been 
used to simplify the implementation of several logical operations for quantum information processing 
[10,12,17,21,30-33]. 

A. Preparation of pseudopure states 

Pseudopure states (PPS) are the staring point for implementation of various quantum algorithms 
by NMR. In a PPS, populations of all the states except one are equalized [16,17]. Among the 
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various methods of creating PPS, we have implemented the method of Pair of pseudopure states 
(POPS), originally suggested by Fung et.al [21]. POPS requires two experiments: (i) equihbrium 
and (ii) a selective population inversion. Subtraction of these two population distributions yields 
zero populations for all the levels except the inverted pair, which in turn has equal positive/negative 
populations. This acts as a pair of pseudopure states. It has been demonstrated elsewhere, that for 
implementation of certain algorithms, POPS is sufficient to begin with [35]. Moreover, in an N-qubit 
system, POPS can also act as a (N-l)-qubit sub-system PPS. 

Here we have chosen to apply (tt) pulse on the fluorine transition 61 which connect 1 00000) <-> 
1 10000). After the (tt)^^ pulse, a gradient pulse is applied to kill any unwanted coherence created by 
the imperfection of the r.f. pulse. The final populations are measured by a small angle (7r/20) proton 
and fiuorine pulses, which convert within linear approximation, population differences into observable 
single quantum intensities. The resulting spectrum is given in figure 6(b). Subtraction of equilibrium 
populations (figure 6(a)) from that corresponding to figure 6(b) yields the pair of pseudopure states 
|10000)(10000| - 1 00000) (00000 1 . The resulting spectrum is given in figure 6(c). Note that the 
transitions 2, 18, 21 and 37 are of opposite sign to that of 61, as they are progressively connected 
to 61 (see figure 5) and transitions 1 and 9 are of the same sign since that they are regressively 
connected to 61. The spectrum obtained after POPS contains only transitions originating from the 
two levels having non-zero populations. This effectively gives the various transitions connected to 
these levels. The same information is also available in the cross-sections of HET-Z-COSY spectrum 
at the frequency of transition connecting the two levels. A cross-section parallel to loi at the lo\- 
frequency of transition 61 from spectrum of figure 2 is given in figure 6(c'). A good match of relative 
intensities of various transitions within each spectrum of 6(c) and 6(c') confirms the creation of 
POPS in 6(c). It may be noted that, while 6(c) is achieved with a selective inversion of transition 
61, 6(c') is a cross-section of a 2D spectrum obtained using non-selective low angle pulses assuming 
perturbation and observation in the linear regime. 

Another method of creating sub-system pseudopure state is; Spatially Averaged Logical Labeling 
Technique (SALLT) [30]. A hard (7r/2) pulse on the protons equahze the populations within each 
subsystem [A and B], with the value differing between the two subsystems. A gradient pulses is 
subsequently applied to kill the coherences created in the process. A transition selective pulse (tt) 
on one of the ffuorine transitions (we have chosen transition 61) then creates subsystem pseudopure 
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states. The populations corresponding to this state are measured by a small angle (7r/20) pulse and 
the resulting spectrum, is given in figure 6(d). The proton spectra of figure 6(c) and 6(d) are nearly 
identical, confirming the creation of subsystem PPS. The populations of POPS and SALLT differ in 
the sense that POPS yields zero populations of all levels except one in each subsystem, while SALLT 
has equal population of all levels (not necessarily zero) except one in each subsystem. The fiuorine 
spectrum in figure 6(d) is therefore identical to the spectrum in figure 6(b). 

B. Implementation of controlled-NOT and controlled-SWAP gates 

Controllcd-NOT gates are essential parts of any computation. A C^-NOT gate will change the state 
of fifth qubit when the other four qubits are in the state 1 1) . Proton transition 60 connects the states 
|11110) and |11111) (see fig. 5). Hence a (7r)^° pulse performs C^-NOT operation. We start from 
equilibrium state and implement the C^-NOT gate. A subsequent gradient pulse kills any coherences 
created by imperfection of the r.f pulse. The output of the gate is stored in the final populations. 
The final populations, measured by small-angle (vr/20) proton and fluorine pulses after implementing 
C^-NOT gate corresponding to |11111) |11110), yield the spectrum given in figure 7(b). Another 
POPS can be created by subtraction of this population distribution and the equilibrium populations. 
Subtraction of figure 7(a) from 7(b) yields the POPS of |11111)(11111| - |11110)(11110|. The 
corresponding spectrum is given in figure 7(c). Creation of the POPS can also be confirmed by 
cross-sections parallel to LO2 of figure 2 at the a;i-frequency of the transition 60, given in 7(c'). 

Controlled SWAP gates arc becoming increasingly popular in quantum information processing 
[36,37]. A C^-SWAP gate swaps the state of last two qubits when the first three qubits are in the 
state Essentially, this gate interchanges the amplitudes between |11110) and |11101). A sequence 
of three pulses {j^f^-iT^f^-ij^f^ achieves this gate. We start from the POPS of 1 11111) (11111| - 
1 11 101) (11101 1 prepared by subtraction of the equilibrium populations from the populations obtained 
by inversion of the proton transition 55 (figure 7(d)). The creation of this POPS is confirmed by the 
cross-section of figure 2 at the frequency of 55, given in figure (d'). The subsequent implementation 
of C^-SWAP gate [(7r)^^-(7r)52.(^)58] interchanges the populations between |11110) and |11101). A 
gradient pulse is applied after each (7r)-pulse to destroy any coherences created by imperfection of 
pulses. The final state is the POPS of |11111)(11111| - |11110)(11110|, and the resulting proton 
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and fluorinne spectra measured by a small angle (7r/20) pulses is given in figure 7(e). A good match 
between 7(e) and 7(c) confirms the implementation of C^-SWAP gate. 

C. Entanglement creation and transfer 

Entanglement is a remarkable property of quantum systems, which has no classical analogue. 
Entanglement is also a key requirement for many quantum information processing protocols. Many of 
these protocols require transfer of entanglement between qubits. Entanglement transfer has become 
specially useful in the cases where one needs to transfer data to the quantum memory (qubits with 
less decoherence that are used for storing data) from the fast processor (qubits have decoherence 
but can perform universal logic operations) in a manner similar to the protocol proposed by Chuang 
et.al. [38]. Entanglement transfer has been demonstrated previously by Cory et.al. in liquid state 
NMR of weakly coupled 4-qubit spin- 1/2 system [39]. Here, we demonstrate entanglement transfer 
in this strongly coupled system by using the last 4-qubits of the 5-qubit system with the first qubit 
acting as an ancillary bit. We start by preparing the 4-qubit subsystem pseudopurc state by SALLT 
(figure 6(d)). We chose to perform further operations in only one of the subsystems (where 1st qubit 
is in state |1)) [domain B]. Hence, in the next paragraph, only the states of the last four qubits are 
given, but they are to be understood as the states with first qubit being |1). 

Starting from the |0000) pseudopure state, uniform superposition of 2nd qubit was created by a 
(7r/2)2 pulse: |0000) |0000) + |0100) (fig 5(a), B domain). The states were tomographed using 
two-dimensional Fourier transform technique [34]. 2D spectra for measurement of all off-diagonal 
elements of the density matrix are given in figure 8. Figure 8(a) contains the totality of the multiple 
quantum (MQ) spectrum of the four protons obtained by a {tt/2) — r — {tt/2) pulse sequence starting 
from equilibrium population distribution. Figure 8(b) on the other hand contains the MQ spectrum 
starting from the state of |0000) -|- jOlOO) and shows the presence of a single quantum coherence of 
this state at frequency of transition '2'( u^^^) along loi. A controlled not gate between 2nd and 3rd 
qubit is implemented by a (vr)^^ pulse. This gate changes the state of third qubit when the second 
qubit is in state |1) and creates an entangled EPR (Einstein-Podolsky- Rosen [41]) pair between the 
2nd and 3rd qubit given by: (]0000) -|- ]0110)). The presence of the double quantum coherence with 
frequency o;^^^ -|- a;^^^^ along loi in the 2D spectrum for tomography of all off-diagonal elements in 
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figure 8(c) and the absence of all other coherences, indicates the creation of the EPR coherence. 
The tomographed density matrix of the B domain, given in figure 8(d), confirm creation of the EPR 
state. A sequence of three transition selective pulses (tt)^^ — (tt)^^ — (tt)^''' swaps the states of |0110) 
and 1 1001), thereby transferring the entanglement to the first and fourth qubit: (|0000) + |1001)). 
Hence the entanglement initially created between the second and third qubit is now transfered to 
the first and fourth qubit, while the second and third qubit are in the |00) state and can be used for 
further computation. All the pulses and the receiver were phase cycled by a standard CYCLOPS 
phase cycle [23]. The presence of the double quantum coherence with frequency lo^'^^ + co^^^^ along 
cui in the 2D spectrum in figure 8(e) and the resulting tomograph of the B domain, given in figure 
8(f), confirm that the entanglement has been transferred with 85% fidelity. 

IV. CONCLUSION 

Residual dipolar couplings can be used to increase the number of qubits for quantum information 
processing. We have used a molecule containing 5 spins oriented in a liquid crystal matrix, and 
exploited the residual dipolar coupling to demonstrate it as a 5-qubit quantum computer. Transition 
assignment is performed using HET-Z-COSY experiment and qubit addressability is achieved by 
transition selective pulses. It is expected that using this protocol higher qubits can be achieved. 
Preliminary experiments promise a 8-qubit system [40]. 

For the implementations reported in this paper, evolution under the internal Hamiltonian was 
not explored. It is however interesting to investigate how effectively the evolution under internal 
Hamiltonian can be manipulated to implement quantum algorithms in these systems. Efforts are 
ongoing in this direction in our laboratory, and recently a new method of preparing pseudopure states 
in oriented systems by exciting selected multiple quantum using evolution under effective dipolar 
Hamiltonian, has been reported [42]. 
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Figure Captions 

Figure 1: Equilibrium spectrum of l-chloro,2-fiouro benzene dissolved in liquid crystal solvent 
ZLI-1132. (a) The spectrum of four strongly coupled protons . (b) The fluorine spectrum. There are 
82 transitions with significant intensities of which 1-60 are proton transitions and 61-82 are fiuorine 
transitions. One fiuorine transition did not show any connectivity in the subsequent HET-Z-COSY 
spectra and could not be assigned to any energy level in the energy-level diagram. This transition 
is shown by asterisk (*). 

Figure 2: Phase-sensitive HET-Z-COSY spectrum of oriented l-chloro,2-fiouro benzene, (a) and 
(c) is obtained when protons are observed in t2 while (b) and (d) is obtained when fiuorine is observed 
in t2. 2Kx2K data points were collected with 8 scans per ti increment. The data was zero-filled to 
4Kx4K before Fourier transform. The spectrum in (a) has the diagonal and connectivity informa- 
tion between proton transitions, while that in (c) has the connectivity information of the fluorine 
transitions to proton transitions, (b) has the diagonal and connectivity information between fiuorine 
transitions, and (d) has the connectivity information of the proton transitions to fiuorine transitions. 
Artifacts due to incomplete suppression of axial and double quantum diagonal are also visible in the 
spectrum . These artifacts are removed on symmetrization (fig. 3) 

Figure 3: Symmetric 2D spectrum of figure 2. A peak picking was performed on the 2D spectra 
of figure 2. The data from the two experiments were added to create a single 2D spectrum. This 
2D spectrum was then symmetrized by a program written in Matlab [28]. The symmetrization is 
performed by retaining the lower absolute value of {uji^ujj) and {ujj^uji) cross-peaks, at both the 
cross-peak positions. For better visibility, the symmetric 2D spectrum is plotted with identical ab- 
solute intensity for all the peaks. The positive peaks are shown by filled circles while the negative 
peaks are shown by empty circles. 

Figure 4: Energy-level diagram of the 5 spin system of oriented l-chloro,2-fiouro benzene. There 
are 32 (2^) energy levels which can be divided into two separate domains, A and B of 16 energy 
levels each. Domain A contains energy levels with fiuorine in -1-1/2 state and domain B contains 
energy levels with fluorine in -1/2 state. Proton transitions connect energy levels within the two 
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domains whereas fluorine transitions connect energy levels between A and B. 

Figure 5: Labeling of the eigenstates of 5-spin system to 5-qubit labels. The first qubit is in state 
|0) in domain A, whereas it is in state |1) in domain B. We have followed a conventional labeling 
scheme here such that in the limit of weak couphng, the state a{+l/2) is given label '0' and the 
state (3{—l/2) is given the label '1'. 

Figure 6: Preparation of pseudopure states in the 5-qubit system of oriented l-chloro,2-flouro ben- 
zene. All spectra were obtained by using a readout small-angle (7r/20) proton or fluorine pulse, (a) 
Equilibrium and ^^F spectra, (b) Spectra corresponding to inversion of population between the 
states |00000) and 1 10000) implemented by a selective (vr) on transition 61. (c) Spectra obtained by 
subtraction of (a) from (b) yielding P0PS(61) of 1 10000) (10000| - |00000)(00000| . Note that the 
transitions 2, 18, 21 and 37 are of opposite sign to that of 61, confirming that they are progressively 
connected to 61 (see figure 5). The transitions 1 and 9 are of the same sign implying that they are 
regressively connected to 61. Transition 19 is of very small intensity and not visible in (d). (c') 
Cross-section parallel to UJ2 of the spectrum in figure 2 at the a;i-frequency of transition 61. (d) 
Spectra corresponding to the 4-qubit subsystem pseudopure state by SALLT (see text for details) . 

Figure 7: Implementation of logic operations in 5-qubit system of oriented l-chloro,2-fiouro ben- 
zene, (a) Equilibrium and ^^F spectra, (b) Spectra corresponding to C^-NOT gate, implemented 
by a selective (tt) on transition 60, which connects jllllO) <-> |11111) (see figure 5). (c) Spectra 
corresponding to the creation of POPS(60) |11111)(11111| - |11110)(11110| by subtraction of (a) 
from (b). In these spectra, transitions 14, 38, 52 and 80 are of opposite sign to that of 60, signifying 
that they are progressively connected to 60, while 55 and 82 are of the same sign confirming that 
they are regressively connected. The transition 45 is of very low intensity and not visible in (c) and 
(c'). 

(c') Cross-section parallel to 102 of figure 2 at the a;i-frequency of transition 60. (d) POPS(55) of 
|11111)(11111| — |11101)11101(| created by a selective (tt) pulse on 55, and a subsequent subtraction 
of equilibrium. Transitions 33, 58,67, 73 and 75 are of opposite sign to that of 55, signifying that 
they are progressively connected to 55 (see figure 5). The transitions 60 and 82 are of the same sign 
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confirming tiiat they are regressively connected to 55. Transition 45 is of very low intensity and not 
visible in (f). (d') Horizontal cross-section of figure 2 at the frequency of transition 55. (e) Implemen- 
tation of C^-SWAP gate on the POPS of (a). A sequence of three selective pulses (7r)^^-(7r)^^-(7r)^^ 
implements the C^-SWAP gate thereby creating the POPS of |11111)(11111| - |11110)(11110|. A 
good match between (e) and (c) confirms implementation of C^-SWAP gate. 

Figure 8: Entanglement transfer within the last 4-qubits of the 5-qubit system of oriented 1- 
chloro,2-flouro benzene. In all 2D experiments a date matrix of IKxlK were collected all zero- filled 
to 2Kx2K before Fourier transform, (a) Starting from equlibrium, multiple quantum spectrum of 
the four protons obtained by a pulse sequence of (7r/2) — r — (7r/2) where r = 30ms, followed 
by the pulse sequence of tomography [34]. (b) Spectrum corresponding to measurement of all off- 
diagonal elements of the state |0000) + |0100). A single coherence along uji at the frequency of 
transition 2,( u^'^^) and absence of any other coherences confirm creation of this state, (c) Spectrum 
corresponding to measurement of all off-diagonal elements of the entangled state |0000) -|- |0110). A 
double quantum coherence along ui at the frequency u^"^^ +u^'^'^^ and absence of any other coherences 
confirm creation of EPR coherences, (d) Corresponding tomographed density matrix using the 
method outlined in [34]. (e) Spectrum corresponding to measurement of all off-diagonal elements 
of the state after entanglement transfer. A double-quantum coherence is observed along uji at the 
frequency u'^'^^ +cj*^^°^ with a small coherence at cu'^'^^ +uj^'^'^\ (f) Corresponding tomographed density 
matrix confirm transfer of entanglement with 85% fidelity. 
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